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We show that as the number of vortices in a three dimen- 
sional Bose-Einstein Condensate increases, the system reaches 
a "quantum Hall" regime where the density profile is a Gaus- 
sian in the xy-plane and an inverted parabolic profile along z. 
The angular momentum of the system increases as the vor- 
tex lattice shrinks. However, Coriolis force prevents the unit 
cell of the vortex lattice from shrinking beyond a minimum 
size. Although the recent MIT experiment is not exactly in 
the quantum Hall regime, it is close enough for the present 
results to be used as a guide. The quantum Hall regime can 
be easily reached by moderate changes of the current experi- 
mental parameters. 

In a recent experiment , the MIT group led by Wolf- 
gang Ketterle has generated as many as 160 vortices in 
a Bose-Einstein condensate of Na atoms. While the vor- 
tices form a hexagonal lattice as expected for equilibrium 
systems, many observed properties appear to differ sig- 
nificantly from the predicted ones 0]. It is thought that 
the system may not be in equilibrium. 

Equilibrium or not, it is clear from the recent Paris 
and MIT experiments [|l| that one can deposit a large 
amount of angular momentum to the condensate by driv- 
ing it with an asymmetric potential at a rotation fre- 
quency fi^ close to the quadrupolar resonance. In view 
of the recent rapid developments ||l| ||], it is conceiv- 
able that higher, or even much higher angular momen- 
tum states can be achieved in the near future. The high 
angular momentum states are interesting because they 
resemble some of the novel solid state systems such as 
type-II superconductor and quantum Hall liquid. How- 
ever, the unique features of Bose-Einstein condensates 
also introduce important differences. 

The purpose of this paper is discuss the properties of 
Bose-Einstein condensates with large number of vortices. 
We show that the physics in this regime is closely related 
to two dimensional quantum Hall physics, even though 
the system is three dimensional. The specific wavefunc- 
tion in this regime allows one to calculate many proper- 
ties of the system, including the effective rotational fre- 
quency = dE/dL. It is clear that fl will not be the 
same as the driving frequency Q^, since it depends only 
on the total angular momentum L, which can be var- 
ied by varying the duration of the drive, even though 
is held fixed. The fact that the vortex lattice has the 
equilibrium form strongly suggests that the system is in 
quasi-equilibrium, characterized by an L and Q which 
change slowly because of the residual asymmetry in the 
confining potential. It is conceivable that some of the 
discrepancy with theory mentioned in ref. Q might dis- 
appear once the effective frequency ft is identified. 

We shall see that in the quantum Hall regime, Corio- 



lis force prevents the unit cell {v) of the vortices lattice 
from shrinking below v* = 7ra^ where a_^_ is the oscilla- 
tor length in the xy-plane. For a large condensate, the 
density profile along z is an inverted parabola similar to 
the stationary case. However, it is a Gaussian in the xy- 

3/2 

plane, with a width a that scales as Rz , where i?^ 
is the maximum extent of the condensate along z. The 
angular momentum of the system is proportional to ttct^ , 
which grows as the size of unit cell v shrinks. 

The quantum Hall regime and the minimum 
vortex separation; The condensate wavefunction 5* is 
determined by the Gross-Pitaevskii functional 
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where h„ =: 
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iMw2^2 and h 
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^Mui'^r'^ are the single particle Hamiltonians along z 
and in the xy-plane, u>^ and a;_L are the frequencies 
of the harmonic potentials in z and in the a;y-plane; 
r = [x, y), = —ihz • r x V is the angular momentum; 
g — Anh ttg^/M, a^^ is the s-wave scattering length, and 
(ri, N) are Lagrange multipliers specifying the angular 
momentum and particle number of the system. It is use- 
ful to rewrite ft.^ — ilL^ = Hl + h{LLij_ — fl)L^, where 



Hl = i-ihV^ - Muj^i X r)^ /(2Af) 



(2) 



Eq.(||) is precisely the Hamiltonian 2It(^*^^ ^ eA/c)^ 
of a charge — e particle moving in xy-plane subjected to 
a magnetic field Bz with a vector potential A — \Bi x r, 
eB/Mc = 2Q,. The eigenfunctions of Ti,^ are 



.\u\'^ / (2a\') Qm Qn ^-\u\^ / a\ 
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with eigenvalues e„ ,„ — huj^[2n + 1) where n and m 

are integers 0, 1, 2, .., u = {x + iy)la^, aj_ — yJKjMu:^, 
d± = {a^/2){dx^idy). The integer n is referred to as the 
Landau level index, and m labels the degenerate states 
within a Landau level n. Since u„ „ are also eigenstates 
of = h{ud^ — with eigenvalue fi{m — n). The 

eigenvalue of — fii^ is h[{Lo^ + -f — VL)m + 
LOj^ . The degeneracy of the Landau level is lifted by the 
deviation uj^ — Vl. [|| 

Interaction effects will mix different (n, m) states. For 
sufficiently small g and uj^ — fJ, the system primar- 
ily resides in the lowest Landau level (LLL), (n = 0), 
since other levels are separated from it by an energy gap 
> At first sight, this might appear to be difficult 

to achieve for it requires gn << hujj_, which implies ei- 
ther g is exceedingly small or density the n is very low 
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Q. While this is true for low angular momentum states, 
it is not true in the high angular momentum limit. As 
we shall see, as the number of vortices increases, (corre- 
sponding to f2 very close to ujj_), the density is thinned 
out sufficiently in the xy-plane so that the LLL regime 
can be attained. We shall therefore proceed by assuming 
that the system in the LLL, and return later to show that 
mixing of higher Landau level is indeed small in the fast 
rotating limit for large condensates with typical interac- 
tions. 

Before proceeding, we note that since the dynamics 
along z is identical to the non-rotating case, we can apply 
the usual Thomas-Fermi approximation along z to ignore 
the kinetic energy jV^^P term and justify it at the end. 
Our task is then to minimize K = J dzlC{z), 

JC{z) = J dr^* [h^ - QL, - /i(z)] ^ + lgj (4) 

with the constraint N = J where fj.{z) = fi — 

^Mlu^z"^. Since iio,m(r) c>c u"^e^^' ^^""^ , a general wave- 
function function in the lowest Landau level is 0(r) = 
f{u)e~'^ /2ai^ where f{u) is an analytic function of u. 
Because of the Guassian cutoff, it is sufficient to take / 
as a polynomial of power Q, which can be written in a 
factorized form na=i (^~^q) according to the fundamen- 
tal theorem of algebra, where are the zeros. Since / 
undergoes a 2-k phase change as u encircles , the factor 
{u — b^) describes a vortex at b^, = (x -I- iy) h^. If 
Q extends to infinite, / will be an infinite product. The 
most general form of "if within the lowest Landau level is 
then •9{x,y,z) = C{z)lUu - b^{z))e-''''^^'^''^-L\ or 



^{x,y,z) = f{z)(f){x,y;z), 



dr|(/)p 



where (j){x, y; z) = y; z)/D{z), 



(5) 
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and = J dr|(/)p is the normalization constant for (j). 
The number constraint now becomes /d.z|/p = N. Not- 
ing that within LLL, 



and eq.(E|) becomes 



(7) 



/C(z) = + h{cu^ - n){r')^) f + (8) 

where /2(z) ~ ^—hfl— ^Mw^z^, {r"^)^ — J r^|(/)pdr and 
I^j, — /dr|0|^. The minimization is performed by varying 
To evaluate (r^) and we note that 

|^|2 = e-«, 7^ = rV«'-2^1n|r-bJ. (9) 



The quantity H in eq.(^ is precisely the energy of a 
charge q = —1 particle in two dimension interacting with 
a uniform positive charge background = tt/o^, and a 
set oi q = —1 negative charges located at {b^}. From 
eq.(||), we have 



Vln - Att ^(Tra^)-! - ^ 5(r - b„) j 



(10) 



which is the Gauss's law for the 2D charged system. 
If {h^} forms an infinite regular lattice, b„ „ = 

n^Ci -I- n2C2, (ni,n2 integers), then J^a^i'^ ~ ^a) = 
u"^ e'^ '', where v = \ci x Cjl is the size of the 
unit cell, and K = ^^Kj^ + £2^2 (^1,^2 integers) are the 
reciprocal lattice vectors; (Ki = (27r/ti)c2 x z, K2 = 
(27r/t))z X Ci). Eq.(^ and (^ then become 

1 1 TT , , 

:3 = -T--. (11) 
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The Gaussian in eq.(|T^) corresponds to replacing the 
vortex density ^aS{r — b^) in eq.(lO) by its average 
1/v. In the 2D electrostatics analog, it corresponds to 
reducing the uniform positive background by the aver- 
age density of the discrete negative charges. Since the 
renormalized background is less confining for the nega- 
tive unit charge at u, we have a > aj_. Eq.(|ll|) also 
requires v > v* — 7ra^, for otherwise 4> will not be nor- 
malized. For hexagonal lattice, this means that the dis- 
tance c between neighboring vortices can not be shorter 
than (27r/A/3)i/2aj_. 

The angular dependence of is contained in the 
product in eq.(|l2|). Since vK.'^ does not depend on the 
length of the basis vector c^, but only on the lattice type 
(i.e. square or hexagonal), the product in eq.(p^) deter- 
mines mainly the lattice structure, and is less important 
in determining the overall density profile of the system. 
We shall therefore proceed by replacing the vortex den- 
sity in eq.(|lC|) by its average, (referred to as "averaged 
vortex approximation" , and come back later to consider 
the effect of the product in eq. ( p^ ) . We shall sec that the 
results of the averaged vortex approximation are intact 
except for changes of numerical factors (of order 1). 

Averaged vortex approximation: Keeping only 
the Guassian in eq.(|l2|), the normalized function in 
eq. 



becomes 4>{r) — -^e ^ . Eq.(fe|) then becomes 



/C(z) 
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fl{z) 



{ajar 

The optimal a and / are 
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9?i^a;_|_(w^ — Vl) 
(14) 
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Since c < aj_, the maximum extent of the cloud along 
z satisfying this condition (denoted as R^) is given by 
3?i(w^ — O) = jl{Rz)- Note that within the average vortex 
approximation, f{z) must vanish for \z\ > in order to 
keep \1/ normalized. We then have 



/i(z) = \mu;1{rI -z^) + shiuj^ - n). 

The number condition N = J f^dz then gives 

1/5 



(15) 
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where (1 + ...) = (1 - 3x + 0(a;2) + ..), x = [(^)^/^(l- 
_0_)3/5]/[i3|iv^]2/5^ As we shaU see, x « 1 tor typi- 
cal experimental parameters and can be ignored. From 
eq.(|l^) and (p[), the size of the unit cell is given by 



v{z) 



6(w_L - ^) 



Rl 
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In the fast rotating limit, fl — > u^, eq.(|Tj) is essentially 
the usual relation — 2Vl/{h/M) for estimating the vor- 
tex density but in a more detailed form. (The connec- 
tion is easily seen since n„ — l/w, and 7ra^ — 'KMuj^l'h « 
2MQ,/h. Note also that even though shrinks quadrati- 
cally as z increases (eq. ([l^) ) , v{z) (and hence the spacing 
between vortices) changes very little except near the tips 
of the condensate at z « i-^z- 

A useful quantity is (Trcr^/u)^ = [ttct^/wJ^^Oi the num- 
ber of vortices inside an area ttct^ at the center cross 
section (i.e. z — Q). From eq.(p^) and (^^, we have 



TTCT 

V 



6(w^ — fi) \a1 



5) (18) 



Combining eq.(|18|) and eq.(|16[), we can express the equi- 
librium frequency f2 and R^ in terms of the number of 
vortices in the center cross section. 
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From cq. ^ and (|l|), we find that the 

al/al = {AbN/^){aJa,){u:Ju:,){ajR,f 



2/3 



(19) 



(20) 



(21) 



The quantities f2 and R^ for systems with different par- 
ticle number and vortex numbers are shown in Table I. 
Finally, we note that the total angular momentum eq.(M) 



of the system is {L^) = (^)- The total energy 
is E ^ JC + ^{LJ + fiN. To the leading term in N, we 
have E^^^[^y {l^ ^)"' (^) (i?JaJ^ hence 
E oc iV^/^ 

Mixing of higher Landau levels and the validity of 
Thomas-Fermi approximation (TEA): The validity of 
TFA requires requires (?iV2M)|V^1'|2 << m(^)I*P, 
which is h^/{2MRl) « fi or {R^/oj^)^ » 1. This 
is easily satisfied for the range of parameters listed in 
Table I. Mixing of higher Landau level will be unimpor- 
tant if the interaction energy density huo ^{a^/ a)'^a^^f'^ 
is much less than that of higher Landau level hujj_f'^, 
or r = {a^^/ a^a^.J'^ « 1. From eq.(|lj), we have 
T = jl{z = 0)l{ZhujJ, or F = ^ {Rja^f. The values 

of F for a variety of external parameters are shown in 
Table I. Although the recent MIT experiment [Q is not 
yet in the LLL regime, it is not too far away. It is clear 
from Table I that the LLL regime can be reached within 
the capability of current technology. Note also that as 
O —t Lo^ and the system is driven more toward the LLL 
regime. However, R^ja^ also decreases. As R^ becomes 
very close to a^,, Thomas-Fermi approximation will not 
be valid and a more accurate treatment along z is needed. 
The present treatment is valid for the frequencies f2 such 
that Rz/ciz > I7 which is satisfied for the cases listed. 

Beyond averaged vortex approximation: A full 
calculation of the quantities (r^)^ and in eq.(^) re- 
quires considering the "structural" product in eq. (p^ . 
We shall show that in the limit oj^, the major effect 
of the structural product is to stabilize the hexagonal lat- 
tice over the square lattice. It does not alter the result 
(r^)^ of the averaged vortex approximation, but changes 
/0 from l/(27rcr2) to a/{2Tra'^) with a - 1. Thus, aU re- 
sults of the simple average vortex approximation remain. 

To prove the above statements, we recall that 
g-acos. = E„=o,i,2,..(-l)"4(")e"', where /„(«) are 
the modified Bessel functions. Eq.(p^) then becomes 



(-l)"K/„^(Cj,)e' 



E 

["k] 



(22) 



AK,Ck] = n (-1)""^"k(Ck), PK] = 5]KnK, 

(23) 



where ^ in eq.(|22|) means the K = term is excluded, 
and [tik, CkI denotes the entire sets {rij^} and {CkI- The 
normalized function is then 



101' 



K#0 ' 



^-CkcosK-t 



(24) 



7rcr2 Z„2 [Ck] 

Z,. [Ck] - E AK, CK]e-'^'^f"-l'/^ (25) 
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We then have (r^)^ 
where 7 = Z„2[(ji]^^ 



2 

(7^7 



1 - 



and a — Z, 



0-2/2 



/|</)|4 = «/(27ra2), 
[2Ck] ^.^[Ck]^^ As 



fi w |^^cr2p[nK]^ >> 1 for all P[nK] 7^ 0. Only terms 
in eq.(p5|) with K?t.k = are important. This im- 
plies 7=1, {r'^)cj, — o"^ , and 



a = Z[2CK]/(^[CK])^ 



^[CK]=En(-i)""^"K(CK), 



(26) 



(27) 



where means summing over with such that 

To evaluate a, we rewrite Z as -^[Ck] = 
W^[CK]E'[LjnK^o(-l)""^"K(CK), where W[C^] = 

IIk^o^oCCk), and /„(C) = /n(C)Ao(C)- Let us de- 
fine Co = Ck ' where Kq is the shortest reciprocal lat- 
tice vector. Since Co = •\/3/(27r) and I/tt for hexago- 
nal and square lattice respectively, and since /n(C) — 

Tjk=o ^i\(l+k)\ ' t^^*^ ^""^^ eq.(|^) are fast convergent. 
Regarding Z formally as a power series of Co j we have for 



the hexagonal lattice, ^"[Ck] = 1 + (3Co/2)/3] 
^[Ck] = W[U 



"^Co r. Co o 



and 



(28) 



where /?„ = E,°r=i E^^^ol^? + ^1 + ^344)-^". The Co' 
term in W comes from expanding /^(Ck) in Co to second 
order in Co- The Co term in eq.(p8|) comes from pairs of I-^ 
terms in Z . The factor (3i is the result of summing over 
all such pairs. The Co term in eq.(p8|) comes from triplets 
of vectors of the same length in eq. (I?!). These triplets 
give rise to a negative sign which lowers the energy of 
the system. The factor /Jj is the result of summing all 
these isosceles triplets. Using eq.(p8|) and eq.(p6|), we 
have a = 1 -|- |Co/3i — ■|Co/^2 + ■■■ ^o orders less than 
Co- The result is a = 1.38, (/3i = 1.218, [32 = 1.001, 
Co = 0.276). 

Repeating the same calculation for the the square lat- 
tice. We have, to the lowest order in Co) W^[Ck] = 1 + 
C„^i + .., Z[Ck] = W^[Ck](l + Ko7i + -), where 7i=E?r=i 
TZ^o i^l + il)'^ = 1-506, which gives a = l + 3Coli + - 
to orders less than Co- The result is a — 1.45, since 
Co = 0.318. The square lattice has energy higher than the 
hexagonal one because it has no negative triplets terms. 

Final Remarks: We have shown that as the number of 
vortices increases, a Bose-Einstein condensate will be- 
come quantum Hall like. A natural question is what 
happens if the angular momentum of the system keeps 
increasing. We have seen that decreases as angular 
momentum increases. Eventually, it will reduce to the 
point where only one harmonic state in the z-direction 



will be occupied. In such "ultra-fast" rotating limit, fluc- 
tuations effects in the xy plane will be important and one 
has to go beyond mean field treatment. The discussion 
of the ultra-fast limit will be presented elsewhere. 
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50 


84 


20 


40 


0.64 


6.8 


1.85 
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50 


1000 


10 


40 


0.179 


10.9 


0.198 
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50 


1000 


10 


100 


0.038 


9.3 


0.146 
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84 


20 


40 


0.13 


4.6 


0.86 
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1000 


10 


40 


0.039 


7.4 


0.092 
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1000 


10 


100 


0.008 


6.38 


0.06 



Caption: Row (a) corresponds to the parameters of the 
recent MIT experiment [|l|. We have counted the vor- 
tices as follows. The Gaussian density profile |</>p in 
eq.(|l^) drops by a factor of 100 from the center to a 
radius r = 2a. If there are 160 vortices within a radius 
of 2(7, then -ko^ jv ~ 40. To remain in the lowest Lan- 
dau level (LLL), one needs F << 1. The parameters of 
recent MIT experiment is not yet in the LLL regime. 
However, the latter can be reached by changing the pa- 
rameters as indicated in case (b) to (f). For all cases, the 
quantity x defined after eq.(p^) is ~ 10" 
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